Exact treatment of Ising model on the helical tori 



Tsong-Ming LiawQ 
Computing Centre, Academia Sinica, 11529 Taipei, Taiwan 

Ming-Chang Huanjj] and Yen-Liang Chou 

Department of physics, Chung-Yuan Christian University, Chungli, Taiwan 

Simon C. Lin 

Institute of Physics, Academia Sinica, 11529 Taipei, Taiwan 
Feng- Yin Li 

Department of Chemistry, National Chung-Hsing University, Taichung, Taiwan 

(Dated: February 2, 2008) 

The exact closed forms of the partition functions of 2D Ising model on square lattices with twisted 
boundary conditions are given. The constructions of helical tori are unambiguously related to the 
twisted boundary conditions by virtue of the SL(2, Z) transforms. Numerical analyses reveal that 
the finite size effect is irrelevant to the chirality equipped with each helical boundary condition. 

PACS numbers: 

Since Onsager obtained the exact solution of the two-dimensional (2D) Ising model with cylindrical boundary 
condition (BC) in 19440, the exact treatments of Ising models on different 2D surfaces have been continuously 
attempted. Most recently, Wu and Lu Q have provided analytical treatments for the Ising models with BCs of 
particular class, including Mobius strip, Klein bottle and self-dual BC. The exact study of the model subject to 
BCs is of fundamental importance. First, it represents new challenges for the unsolved lattice-statistical problems 
UEIlllSli0EI3[l3- Second, it is crucial for the finite-size analysis [ill H 13 III E3 Furthermore, it 
provides an optimal testbed for the predictions of the conformal field theory 0] . Numerical simulations are plausible 
for the exact analyses and provide very rich content for the theory of finite-size scalings[T^|. For example, based on 
the exact analysis of dimer statistics, by Wu and Lu (0, 1998), Kaneda and Okabe 01 have achieved, via computer 
simulations, more thorough understanding for the finite-size scaling behaviour of the Ising models subject to the 
boundary types of Mobius strip and Klein bottle. While interesting numerical studies, concerning the excess number 
of percolation^] and the Binder parameter^f|, for the Ising model for the twisted BCs further proceed, the problem 
for lacking the corresponding closed form of the partition functions turns out to significant. 

Boundary conditions are prescribed by sets of primitive vectors which impose the periodicity on the corresponding 
directions. For definite BC, sets of primitive vectors are by no means unique [Tsj ]. For 2D, the equivalent transfor- 
mations among the primitive vector-pairs on lattice essentially preserve the area spanned by the vector-pairs and are 
thus recognised as SL(2, Z). This is the prototype of the modular symmetry discussed in the context of conformal 
field theory [lfj. The helical BC is of particular significance owing to its geometrical feature and its relevance to the 
formulation for the nanotube physics [TJ- Helical tori are formed by pairwise joining the edges of the sheet spanned 
by any two orthogonal primitive vectors. The construction ends up with distinct orientations of the underlying lattice, 
labelled by the chirality as well as the chiral aspect ratio. The conventional periodic BC is referred as the helical 
BC with trivial chirality, as depicted in Fig. 1. The twisted BC, on the other hand, counts on the modification to 
the conventional BC by cutting the torus and then rejoining after twisting. Furthermore, the helical BC is shown to 
be the subclass of twisted one according to the equivalence relations, as depicted in Fig. 2. 

In this Letter, the 2D Ising model subject to the twisted BCs is exactly analysed. The general form of such 
partition functions is obtained, firstly. Symmetry conditions are employed to reduce the redundancy on setting the 
twisting factor a in relation to the conventional aspect ratio A. In addition, any helical BC is shown unambiguously 
equivalent to a definite twisted BC, by virtue of the SL(2, Z) transform. The finite-size shift of critical temperature 
is thus investigated numerically. It turns out that the scaling behaviour is found chirality-independent. Meanwhile, in 
examining the twisting pair dependence, the A — 1 scaling rule appears to be twisting-independent. We then conclude 
by few remarks on the comparison to the previous numerical issues. 

Consider a M x N square lattice with the coordinates of the lattice sites specified in form of xm + yn. The partition 
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function of the Ising model lattice is given as Zm,n = [ 2cosh(/3Ji) cosh(/3J 2 ) ] MN Qm,n with the reduced partition 

function Qm,n = Yl m =iYl n =i ^Qm,n, where Q m ,n = \ H{a mn }\0- + ticr m ,n^m+i,n)(^ + *20'm,nO'm,n+i)] ■ Here we 
use the notations, U = tanh(/3Ji) with Jj, for i = 1,2, denoting the coupling constants along x and y directions, 
and /? = 1/kgT. Twisted i?Cs amount to the identifications of the spin variables whose locations are related by the 
pair of primitive vectors, say, a i and a 2 . Basically, only two types of twisting should be in considerations: One 
is referred as Twi(M, N,d/M) specified by the primitive vectors { a\ = Mx + dy, a 2 = Ny}, and the other is as 
Tw n (M, N, d/N) specified by {~a 1 = Mx, ~a 2 = dx + Ny}. 

According to Plechko ( 8], 1985), the reduced partition function takes the form of 
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with 



q m ,n = e n n 

{a mrl } m—1 n—1 
V&W — I dn dn* p a ™." a m,n[ A A* 1 

= J db m , n db* m<n e b — b ^[B m . n B* mjl+1 } (2) 

where A m n = ^ J ra m ^ n U m ^ n ,A^ nn = 1 "H^l^m— l.n^ro.ra) ^m,n = 1 + b m ^ n t7 m ^ n and = 1 + i2^m,n-l f7 m,nj I n 

above, two pairs of conju gate Grassman variables, {a m ^ n , a* n „} and {6 mj „, 6^ „} have been introduced. As technically 
known to the Refs. |3.MllO|. the handling of the boundary Boltzmann weights, 
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*r = n *s« n *Sv, (3) 
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remains central in the treatments. It turns out to be instructive to reexamine the paradigm which solves this problem 
in the original periodic settings a m+M .n = cr m ,n an d a m ,n+N = v m ,n- 

InRef. y| (1985), the boundary Boltzmann weights are rearranged such that \&r = ^7 ki +^7 ^2+^7^3— ^j\r 4 
subject to the BCs I^s, imposed on the Grassman variables, with 
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*7 = / fl II ~B*rn.l fl A M , n fl ^m.JV, (4) 

n m n m 

where the arrows indicate the ordering for the multiplications and we employ the notation J for all the coming weighted 
integration over relevant Grassman variables. Subsequently, mirror ordering is applied routinely and furnishes the 
simple expression of pure Grassmanian integrations, 

Qm,n = i[Gf|r a +G|r 2 + G|r 3 -G|rJ, (5) 

M,N 

G = / exp{ E k.t.^m.ti + *l*20m_i 1 n&m ) n-l + (*l a m-l,n + *2&m,n-l)( a m,n + &m,n)]}, (6) 
m.n 

where the integrations can be diagonalised and carried out, straightforward @- 

The reviewing paragraph above suggests that the modification is only essential for the twisted BC in the the key 

steps, i.e., from Eq. J2J to Eq. Q). For Twi, the \&^f term in Eq. © is calibrated in relation to the toroidal one. 
This then leads to 



. M > N — d < N < N ^ M < 

*r = / jj jj (1 tioM.fc+do-i,*) n ^ ~ ^^fc+d-jv^i.fc) n n b ^,n, (7) 

J ~_i fc=1 fe=n-d+l — 1 1 



where the -BC Cm+M.n+d = G m.n has been explicitly employed. Reordering of the first three products in Eq. Q 
is essential such that the form of Eq. Q can be achieved. By recursive use of the identity for the permutation of 
Grassmanian functions p|, we employ, instead, 

xyz = ^(zyx- - z-y-x- + zy-x + zy~x-), (8) 
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where X, Y and Z stand for the corresponding three objects and the superscript "— " denotes flipping the sign of 
the Grassman variables. Accordingly, the form of Eq. Q is achieved which implies Eq. (JSJ. Note that the form of 
Eq. @ preserves under twisting. However, the BC's imposed on the Grassman variables are modified in response 
to the corresponding sign flipping appearing in the deduction of Eq.©. For convenience, the compact notation as 
Ti = (±,±) can be employed as follows. The first sign in the parenthesis corresponds to N — ± a* n and the 
second one is for a* M n+d = ± aJ5„. The BCs are given as Ti = (—,—), T2 = (+,—), T3 — (— ,+) and = (+,+). 
The exact partition function is straightforward, henceforth. 
For Twi, with a = d/M, the reduced partition function is 



Qm,N = Ul^M^+^A^+^M^-^^^^M ' ) 
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where Ao = (1 + tf)(l + t^), Ai = 2ii (1 — t\) and A2 = 2t 2 (1 — t\). In addition, the function sgn(x) denotes the sign 
of the value x and T c is the critical temperature of the bulk system. The reduced partition function for Twn remains 
formally as Eq. (J5J) but with 
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where, instead, a = d/N. 

It can be checked for Eq. © that N = Qjf N based on either Eq. (|10|l or Eq. while, intuitively, twisting 

either clockwise or counterclockwise is not classified by the system. Noteworthy is also that reversing the sign of 
a twist factor a can not be obtained via the SL(2, Z) transform. On employing this transform explicitly, pairs of 
primitive vectors are related among each other in the manner of 
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= M 



a 1 
a 2 



V M e SL(2,Z). 



(12) 



Consider Twi, for example. The choice of matrix elements M.\\ = 1, .M12 = J G Z, M21 = and M.22 = 1 gives rise 
to the new pairs of primitive vectors { a 1 = Mx + (d + N)y, a 2 — Ny}, which prescribes the same BC. As evidence, 
Qm n — Q°\i n A can b e explicitly checked, where the conventional aspect ratio appears as A = N/M. Therefore, 
the effective range of a is < a < A. In addition, the equivalence Twi(M, N,a — A/ J) = Twu(M = JM,N — 
N/J,a — 1/a) can be achieved by virtue of choosing the elements A^n = J € Z, M12 = —1, A^2i = 1 and A^22 = 
with a = A/ J and a = J/ A. Again, the partition functions based on Eqs. IjlOl and 111|) appear to fulfil the these 
relations. Hence, it is sufficient to study the unique correspondence of a helical torus to the one of the above twistings, 
say Twj. 

The helical tori, on the hand, lie in the orthogonal primitive vector pair, 



= x Pi 



y Q\ 



= -x Q2 + y P2 



(13) 



where the two radii for the torus are given as Li — y/ Pf + Q\ for i = 1,2. Hence, let the helical system denoted by 
Hl(B, Li, x), where the chiral aspect ratio B = L2/L1 and the chirality \ = Q\/P\ = Q2/P2- In order to furnish 
the equivalent structure Hl(B, L\,x) = Twi(A, M, a), Ain = Pi/M and A^2i = —Q2/M implies that 
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(14) 
(15) 
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However, few remarks on the uniqueness of the relations above remain essential. 
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The [.M12, M22] pair is unambiguously determined up to Ain and M21 for < a < A. This is because shifting 
[7^12,-^22] by appending [J.M11, J.M21] leaves Ea. (fPo)l invariant VJ G Zbut only causes a deviated by J A in Eq. 
p7(l . Meanwhile, the allowable region for [.M12, .M22] appropriate for < a < A is exactly one vector section 
[Alii, .M 21]. Iu addition, the ambiguity relating to size dependence can be removed by noting the coprime properties 
between A4n and .M21 which ensures the solubility of the integer pair [^12,^22] subject to Eq. I|15l) . Consequently, 
a helical two-tuple (B,x) is equipped a unique pair {A, a} for the twisting in the effective range. Moreover, the 
classification of helical via the twisting parameters remains unambiguous. This is because if different helical tori 
were equivalent to the same twisted BC, SL(2,Z) transforms would have been held among them, which can be shown 
impossible. The effective range of the helical BC can also be further reduced. The partition function is unable to 
classify the rolling up direction in forming the tori, hence, no distinction between the characterisations x an d — X 
is essential. Subsequently, once HI (B, \) — Twj(A, a), one can derive that Hl(l/B, l/x) — Twi{A,—a). To be 
concrete, assuming X > 0, x' = l/x implies M' n = —M2i,M' 21 = —M-u as well as M' 12 = M22 and M' 22 = M12 
according to Eq. JT2}. This then preserves A but gives a' = —a by virtue of Eqs. (|15() and i|16|) . 

The shift of the specific-heat peak T max away from the critical temperature T c under the isotropic couplings can 
be computed from the the exact partition function. Upon using the parametrisation in terms of Twi(A, M, a), the 
critical shift 6(A, a) = (T max (a) -T c )/T c is plotted against l/L = I/VAIN in Fig. 3. For helical BCs, Eqs. lTH)) - ltT7|) 
are employed in order to determine the critical shifts 8(B,x) versus l/L = l/\ / L 1 L 2 = 1/VMiV for various Bs and 
XS, as shown in Fig. 4. In Fig. 3, the curves of scaling for definite A deviate by altering the twist factor a. However, no 
such splitting is found for the exceptional case A = 1. On the other hand, all the critical shifts 8 (B, x) with the same 
B value, in Fig. 4, fall into to one single smooth curve, the finite-size effect turns out to be chirality- independent. 
Hence, the matching of the two particular curves, i.e. for B = 1 in Fig. 3 and for A = 1 in Fig. 4, appears to be 
an additional feature. Moreover, the critical shift 8(B,x) flips its sign at B = bo and l/bo with 60 — 3, as it was 
anticipated by Ferdinand and Fisher for the conventional periodic BC, where the exact bo value was determined 
as bo — 3.13927.., a result which now applies for all the helical tori. 

In conclusion, we provide the complete description for the finite-size effect of Ising Model subject to the subclass 
helical BCs of the twisted tori. This is explicitly done by solving the exact form of the partition function appropriate 
for all the twisted BCs. The evidence of the finite-size effect being chirality-independent basically supports the 
invariance of the scaling behaviour of partition function under rotation of the primitive vector pair subject to BCs, 
conjectured in Ref. |15|. However, the particular coincidence for A = 1 and B = 1 regardless a and x suggests further 
interesting points exceeding beyond the rotational invariance. For consistency, we stress the fact that A = 1 does not 
non-trivially permits any helical structure, as one may observe in Eqs. (|14[l - i|17[) . As the final remark, the invariant 
aspect ratio A/(l + a 2 ) 4,9] coincides with B only for x — a , nor does it pertain to the case where A = 1. 

This work was partially supported by the National Science Council of Republic of China (Taiwan) under the Grant 
No. NSC 93-2212-M-033-005. 
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FIG. 1. The formation of helical tori by pairwise joining the edges of the rectangle spanned by any orthogonal set 
of vectors on the lattice plane: (a) the direction of the primitive vectors coincides with the lattice orientations for the 
conventional toroidal BC and (b) the helical tori are formed for the non-coincidence. 

FIG. 2. Equivalence between the BCs in helical and twisted schemes prescribed by {a\, 0,2} and {a[,a' 2 } respectively, 
on a M x N square lattice. For the helical BC, the setting Q1/P1 = Q2/P2 ensures that the two primitive vectors 
are orthogonal. On the other hand, twisting is generated by a (i-unit traverse shift. 

FIG. 3. Plotting 9(A,a) against l/L for A = 1,2,3,4 with a = 0(o), 0.1.4(A), 0.2A(v), 0.3A(o), 0AA(») and 
0.5 A(x). The scaling behaviours are obviously deviated by a. Nonetheless, for A = 1 no splitting is found with 
respect to the twisting factors. 

FIG. 4. The plot of 9(B, x) versus l/L. For a given chiral aspect ratio B. Results of different chiralities x collapse 
into one curve and the curves of both 9(B,x) and 9(1/ B, 1/x) versus l/L coincide. 
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